The balance between excitation and inhibition controls the temporal organization of 

neuronal avalanches 
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Neuronal avalanches, measured in vitro and in vivo, exhibit a robust critical behaviour. Their 
temporal organization hides the presence of correlations. Here we present experimental measure- 
ments of the waiting time distribution between successive avalanches in the rat cortex in vitro. This 
exhibits a non-monotonic behaviour, not usually found in other natural processes. Numerical sim- 
ulations provide evidence that this behaviour is a consequence of the alternation between states of 
high and low activity, named up and down states, leading to a balance between excitation and inhi- 
bition controlled by a single parameter. During these periods both the single neuron state and the 
network excitability level, keeping memory of past activity, are tuned by homeostatic mechanisms. 

PACS numbers: 05.65.-l-b, 05.45.Tp, 89.75.-k, 87.19. L- 



Spontaneous neuronal activity can exhibit slow oscil- 
lations between bursty periods, or up-states, followed by 
substantially quiet periods. Bursts can last from a few to 
several hundreds of milliseconds and, if analysed at a finer 
temporal scale, have often shown a complex structure in 
terms of neuronal avalanches. In vitro experiments record 
avalanche activity from mature organotypic cultures 
of rat somatosensory cortex where they spontaneously 
emerge in superficial layers. The size and duration of neu- 
ronal avalanches follow power law distributions with sta- 
ble exponents, which is a typical feature of a system in a 
critical state, where large fluctuations are present and the 
response does not have a characteristic size. The same 
critical behaviour has been measured in vivo from rat cor- 
tical layers during early post-natal development^, from 
the cortex of awake adult rhesus monkeys [4|, using micro- 
electrode array recordings, as well as for dissociated neu- 
rons from rat hippocampus]^, Q or leech ganglia^]. In 
vitro, quiet periods measured between bursts, also called 
down-states, can last up to several seconds. The emer- 
gence of these down-states can be attributed to various 
mechanisms: a decrease in the neurotransmitter released, 
either due to the exhaustion of available synaptic vesicles 
or to the increase of a factor inhibiting the release such 
as the nucleoside adenosine |8|, the blockade of receptor 
channels by the presence of external magnesium^, or 
else spike adaptation[l3|. A down-state is then character- 
ized by a disfacilitation, i.e. absence of synaptic activity, 
of a large number of neurons causing long-lasting returns 
to resting potentials ll| . Recently, it was shown analyti- 
cally and numerically that critical behaviour 12| charac- 
terizes up-states, whereas down-states are subcritical 13 1. 



Whereas action potentials are rare during down-states, 
small amplitude depolarizing potentials, reminiscent of 
miniature potentials from spontaneous synaptic release, 
occur at higher frequencies. The non-linear amplification 



of small amplitude signals contributes to the generation 
of larger depolarizing events bringing the system back 
into the up-state as observed in cortical slabs dis- 
sociated cultures [ist and slice cultures The analysis 
of the amount of time striatal spiny neurons [ItI, [i3 and 
cortical pyramidal neurons [3| spend at each value of the 
membrane potential shows that both cell types toggle be- 
tween two preferred values [l^: a very negative one in the 
down state, and a more positive, depolarized one, in the 
up-state. The up-state being just a few millivolts from 
the action potential threshold, suggests that during the 
up-state neurons respond faster and more selectively to 
synaptic inputs. For cortical neurons the up-state would 
be a metastable state, i.e. the membrane potential would 
soon decay down to the resting potential value, if network 
mechanisms would not sustain the activity. The up-state 
has therefore network, rather than cellular, properties. 

Here we focus on the temporal organization of neu- 
ronal avalanches both in organotypic cultures and neu- 
ronal networks simulations. Each avalanche i is charac- 
terized by its starting and ending times, and t{ . The 
temporal organization is analysed by evaluating the dis- 
tribution of waiting times Ati — t\^^ — t{ . This is a fun- 
damental property of stochastic processes, widely investi- 
gated for natural phenomena and able to discriminate be- 
tween a simple Poisson and a correlated process. Indeed, 
in the first case the distribution is exponential, whereas 
it exhibits a more complex behaviour with a power law 
regime if correlations are present. For a wide variety of 
phenomena, e.g. earthquakes and solar flares [2l[, human 
dynamics [2§|, biological systems H^, etc., this distribu- 
tion always shows a monotonic behaviour. Recent results 
on freely behaving rats provide a lognormal size distribu- 
tion and a monotonic waiting time distribution uniquely 
controlled by the average occurrence rate 2J| . Anaes- 
thetized rats, conversely, exhibit a more heavy-tailed size 
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FIG. 1: (Color online) The distribution of waiting times for 
seven different slices of rat cortex. All curves show an initial 
power law regime between 10 and about 200ms, with an aver- 
age exponent 2. 15 ±0.32. For At > 200ms curves can become 
quite different with the common characteristics of a local min- 
imum located at 200ms < Atmin < Is, followed by a more or 
less pronounced maximum at At ~ 1 — 2s. In the insets two 
temporal sequences of neuronal activity for numerical (sum 
of potential variations) and experimental {fjV) data. 



distribution and no universal scaling for the waiting time 
distribution. Here we show that the waiting time distri- 
bution for neuronal avalanches in vitro has an unusual 
non monotonic behavior. Numerical simulations on neu- 
ronal networks suggest that this is controlled by the slow 
alternation of up and down states, which determines both 
the network and the single neuron behavior. 

Experiments were performed on coronal slices from rat 
dorsolateral cortex (postnatal day 0-2; 350 fim thick) at- 
tached to a poly-D-lysine coated 60-microelectrode ar- 
ray (MEA; Multichannelsy stems, Germany) and grown 
at 35.5 °C in normal atmosphere and standard cul- 
ture medium without antibiotics for 4-6 weeks before 
recording. Avalanche activity was measured from cortex- 
striatum-substantia nigra triple cultures or single cor- 
tex cultures as reported previously [l]. Spontaneous 
avalanche activity is recorded outside the incubator in 
standard artificial cerebrospinal fluid (ACSF; laminar 
flow of 1 ml/min) under stationary conditions for up 
to 10 hrs. The spontaneous local field potential (LFP) is 
sampled continuously at 1 kHz at each electrode and low- 
pass filtered at 50 Hz. Negative deflections in the LFP 
(nLFP) are detected by crossing a noise threshold of -3 
SD (SD stands for standard deviation ^ 3 — 5fiV) fol- 
lowed by negative peak detection within 20 ms. nLFP 
times and nLFP amplitudes are measured. Neuronal 
avalanches are defined as spatio-temporal clusters of 
nLFPs on the MEA[2^. A neuronal avalanche consists 
of a consecutive series of time bins of width St that con- 
tain at least one nLFP on any of the electrodes. Each 
avalanche is preceded and ended by at least one time 
bin with no activity. The waiting time At is simply 



given by the number of empty bins between two suc- 
cessive avalanches times St. Without loss of generality, 
the present analysis is done with St estimated for each 
culture as the average inter nLFP interval on the array 
and ranged between 3-6 ms for all cultures. 

In Fig.l we show the waiting time distribution for dif- 
ferent cultures of rat cortex slices. The curves exhibit 
a complex non-monotonic behaviour with common fea- 
tures: an initial power law regime and a local minimum 
followed by a more or less pronounced maximum. The 
presence of a power law implies that avalanche occur- 
rence is not a pure Poisson process, namely successive 
avalanches are temporally correlated [26]. Moreover, the 
non-monotonic behaviour is not usually observed in nat- 
ural phenomena. In order to investigate the origin of 
this behaviour, we simulate avalanche activity by a neu- 
ronal network model 27-2^, which is able to reproduce 
the scaling properties of neuronal avalanches. Here we 
question whether and how the complex temporal orga- 
nization of avalanches can be caused by the slow alter- 
nation between up-states and down-states. The basic 
idea is that after a large avalanche the involved neurons 
become hyperpolarized and the system goes into a down- 
state. Conversely, after a small avalanche active neurons 
remain depolarized, the system stays in an up-state. 

We consider N neurons at random positions, charac- 
terized by their potential Vi . Neurons are connected by a 
classical scale- free network [30| , where a neuron i has an 
out-going connectivity degree kouti ■ Once the network of 
output connections is established, we identify the result- 
ing degree of in-connections, km for each neuron. To each 
synaptic connection we assign an initial random strength 
gij, with gij ^ gji, and to each neuron an excitatory 
or inhibitory character, with 10% inhibitory synapses. 
Whenever at a given time the value of the potential at a 
site i is above a certain threshold, Vi > fmax, the neuron 
sends action potentials which arrive to each of the kouti 
pre-synaptic buttons. As a consequence, the total charge 
entering the connected neurons is qi oc Uikouti, as in a 
firing rate based charge distribution. Each neuron re- 
ceives charge in proportion to the synaptic strength gij, 



where the sum is on all 



v,{t + l)=v,it)±Si-^-^, 

out- going connections of i. Here the membrane poten- 
tial variation is obtained by dividing the received charge 
by the surface of the soma of the post-synaptic neuron, 
proportional to the number of in-going terminals kinj ■ 
The plus or minus sign is for excitatory or inhibitory 
gij, respectively. After firing, a neuron is set in a re- 
fractory state lasting 1 time step (about 10 ms), during 
which it is unable to receive or transmit any charge. At 
the end of an avalanche, we implement Hebbian plas- 
ticity rules: The strength of the used connections be- 
tween active neurons is increased proportionally to the 
activity of the synapse [sij, namely the membrane po- 
tential variation of the post-synaptic neuron, gij{t+ 1) = 
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9ij (0 + ("j + 1 ) ~ ''^j (*) ) / ^max ■ Conversely, the strength 
of ah inactive synapses is reduced by the average strength 
increase per bond, Ag = t^9iA^)/-^bi where Ni, is 
the number of bonds. The presence of both strengthen- 
ing and weakening rules implements a homeostatic reg- 
ulatory mechanism for synaptic strengths, which under- 
lies the system's critical behaviour. An external stimulus 
triggers further activity in the system: at the end of each 
avalanche the potential of a random neuron is increased 
by a small amount until another neuron gets at thresh- 
old and starts an avalanche. We implement the plastic- 
ity rules during a series of stimuli in order to modify the 
synaptic strengths, initially random. Previous studies 
have verified that the critical behaviour of avalanche dis- 
tributions does not depend on parameter values or net- 
work properties and that this model reproduces quan- 
titatively the background spectrum of measured EEG 
signals 27 - 2^ . The implementation of a scale free net- 
work of connections in the present study is motivated by 
numerical convenience in terms of cpu time. 

In order to implement the alternation between up and 
down-states, at the end of each avalanche we measure its 
size in terms of the sum of depolarizations Svi of all active 
neurons, sav — X^'^^i- ^^^^ avalanche is larger 

than a threshold, sav > Sau"' system transitions 
into a down-state and neurons active in the last avalanche 
become hyperpolarized proportionally to their previous 
activity, namely we reset 



hSvi 



(1) 



where h > 0. This rule introduces a short range memory 
at the level of a single neuron and models the local inhi- 
bition experienced by a neuron, due to spike adaptation, 
adenosine accumulation, synaptic vesicle depletion, etc. 

Conversely, if the avalanche just ended has a size sav < 
SaI", the system either will remain, or will transition into 
an up-state. All neurons firing in the previous avalanche 
are set to the depolarized value 



(2) 



The neuron potential depends on the response of the 
whole network via sav, in agreement with measurements 
of the neuronal membrane potential which remains close 
to the firing threshold in the up-state, s^™ controls the 
extension of the up-state and therefore the level of ex- 
citability of the system. The high activity in the up- 
state must be sustained by collective effects in the net- 
work, otherwise the depolarized potentials would soon 
decay to zero, and therefore the random stimulation in 
the up-state has an amplitude that depends on past ac- 
tivity. Eqs. (1) and (2) each depend on a single pa- 
rameter, h and s^™, which introduce a memory effect 
at the level of single neuron activity and the entire sys- 
tem, respectively. In order to reproduce the behavior 
observed experimentally, the parameters s^*" and h are 




FIG. 2: (Color online) Waiting time distributions measured 
experimentally are compared with the average numerical dis- 
tributions for 100 networks with TV — 64000 neurons. Top: 
numerical curve (s^*" — 140 and h = 0.017) fitting the exper- 
imental curve with blue squares in Fig.l; Bottom: numerical 
curve (sa™ = HO and h = 0.02) fitting the experimental 
curve with red diamonds in Fig.l. In the insets the waiting 
time distribution evaluated separately in the up and down- 
state for the numerical (upper insets) and the experimental 
curves (lower insets). For the numerical curves, statistical 
error bars not shown are comparable to the symbol size. 



controlled separately. Our simulations will show that the 
ratio R = h/s^™ is the only relevant quantity controlling 
the temporal organization of avalanches. 

Numerical simulations show that the system indeed 
switches between up and down states, with different tem- 
poral durations (insets of Fig.l). The numerical waiting 
time distributions (Fig. 2) exhibit the non-monotonic be- 
haviour of the experimental curves, where the position of 
the minimum is controlled by the value of Sa™ ^^'^ 
power law regime scales with the same exponent ~ — 2 
as experimental data. The agreement between the nu- 
merical and the experimental distribution is confirmed 
by the Kolmogorov-Smirnov test at a p = 0.05 signifi- 
cance level. Both distributions pass the statistical test 
with p = 0.99 (bottom panel) and p — 0.68 (top panel). 
The different contribution from the two states is reflected 
in the activity temporal scale (insets of Fig. 2). The up- 
state generates strongly clustered avalanches, originating 
the power law regime of the waiting time distribution, 
whose extension depends on Sai>"- Large At between 
avalanches generated in the upstate are observed with a 
very small probability, which increases with decreasing 
h. Conversely, the waiting time distribution evaluated 
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in the down-state has a bell-shaped behaviour centered 
at large intertimes which depends on /i, i.e. for a larger 
disfacilitation of the network the probability to observe 
intermediate waiting times decreases in favour of long At. 

The presence of the minimum and the height of the 
relative maximum are sample dependent (Fig.l) and for 
each sample the agreement between numerical and ex- 
perimental data depends on the subtle balance between 
excitation and inhibition. For different samples, optimal 
agreement is realized when the ratio R = h/s^^J^ ~ lO"**. 
For instance, enhancing excitation, by increasing the 
threshold value s^™, clearly produces a major shift in 
the data (Fig. 3). Increasing inhibition, by increasing the 
parameter h, generates the opposite effect, recovering the 
good agreement with experimental data. Interestingly, 
the avalanche size and duration distributions also repro- 
duce the experimental scaling behaviour for the parame- 
ter values expressing the balance between excitatory and 
inhibitory components. The abrupt transition between 
the up and down-state, controlled by a threshold mech- 
anism, generates the minimum observed experimentally. 
Simulations of up-states and down-states only in terms 
of different external drives, without the single neuron 
state dependent behaviour (Eq.s 1-2), provide a mono- 
tonic waiting time distribution (inset Fig. 3). 

This complex non-monotonic behaviour, controlled by 
the system balance level between excitation and inhibi- 
tion expressed by the parameter R, does not simply de- 
pend on the occurrence rate. The different behaviour 
with respect to alive rats [2J] could be attributed to a 
larger separation in characteristic temporal scales be- 
tween up and down-states. Indeed, long lasting down- 
states in our case originate waiting times one order of 
magnitude longer than for awake rats. Avalanches are 
temporally correlated in the up-state, whereas down- 
states are long term recovery periods where memory of 
past activity is erased. A detailed analysis of power spec- 
tra may shed further light on the temporal features of 
this alternation. The good agreement with experimen- 
tal data indicates that the transition from an up-state 
to a down-state has a high degree of synchronization. 
Moreover it confirms that alternation between up and 
down-states is the expression of a homeostatic regulation 
which, during periods of high activity, is activated to con- 
trol the excitability of the system and avoid pathological 
behaviour. The model suggests that the crucial feature of 
this temporal evolution is the different single neuron be- 
haviour in the two phases: These collective effects must 
be supported by the single neuron behaviour, which tog- 
gles between two preferential states, a depolarized one in 
the up-state and a hyperpolarized one in the down-state. 
The model suggests that the depolarized neuron state is 
a network effect: the avalanche activity itself determines 
how close to the firing threshold a neuron stays in the 
up-state. Conversely, the hyperpolarized state is a form 
of temporal auto-correlation in the neuron activity. The 
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FIG. 3: (Color online) Waiting time distribution measured 
numerically for 100 networks of A'^ = 16000 neurons with 
different s^l" and h. The best agreement is obtained for 
R ~ 10~*. In the inset waiting time distributions obtained 
with different stimulations in the up and down-states and 
without the single neuron state behaviour {h — 0). 



critical state realizes the correct balance between excita- 
tion and inhibition via these self-regulating mechanisms. 
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